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Ch. Hartinger, F. Mayr, and A. Loidl 
EP V. Center for Electronic Correlations and Magnetism, University of Augsburg, 86135 Augsburg, Germany 



in 

o 
o 

(N 

CD 

m 






I 

o 



> 

m 

o 
o 



o 
o 



- '—I 

X 



T. Kopp 

EP VI, Center for Electronic Correlations and Magnetism, University of Augsburg, 86135 Augsburg, Germany 

(Dated: February 2, 2008) 

In the colossal magnetoresistance manganites polarons have been proposed as the charge carrier state which 
localizes across the metal-insulator transition. The character of the polarons is still under debate. We present 
an assessment of measurements which identify polarons in the metallic state of La2/3Sri/3Mn03 (LSMO) 
and La2/3Cai/3Mn03 (LCMO) thin films. We focus on optical spectroscopy in these films which displays a 
pronounced resonance in the mid-infrared. The temperature dependent resonance has been previously assigned 
to polaron excitations. These polaronic resonances are qualitatively distinct in LSMO and LCMO and we 
discuss large and small polaron scenarios which have been proposed so far. There is evidence for a large 
polaron excitation in LSMO and small polarons in LCMO. These scenarios are examined with respect to further 
experimental probes, specifically charge carrier mobility (Hall-effect measurements) and high-temperature dc- 
resistivity. 

PACS numbers: 75.47.Lx, 71.38.-k, 72.80.-r, 78.20.-e 



I. INTRODUCTION 

Since the discovery of the colossal magnetoresistance ef- 
fect in thin films of La2/3Cai/3Mn03 by Helmolt et al.,^ the 
properties of doped perovskite manganites have been a focal 
point of research in correlated electronic systems. These com- 
pounds reveal a variety of ordinary and exotic ground states, 
for example spin, orbital or charge order, cooperative Jahn- 
Teller (JT) distortions or electronic phase separation, as a re- 
sult of a delicate interplay of orbital, spin, charge, and lat- 
tice degrees of freedom. While in the early days research 
concentrated on the appearance of colossal magnetoresistance 
(CMR), more recently also issues like orbital order, the impor- 
tance of electron-phonon coupling, or phase separation sce- 
narios promoted among the topics of interest. Optical spec- 
troscopy, in both single crystals and thin films, has contributed 
to unravel the complex physics and provided insight into the 
importance of electron-phonon coupling for modelling the op- 
tical conductivity in the manganites. 

Okimoto et aL^ and Kaplan et al? reported on a large 
transfer of spectral weight from high to low energy with de- 
creasing temperature for Lai_a;Sra;Mn03 single crystals and 
for Ndo.7Sro.3Mn03 thin films, respectively. The observed 
temperature dependence is consistent with the theoretical ap- 
proach by Millis et al.* who went beyond Zener's double 
exchange^ and included dynamic JT polaron effects. Sub- 
sequent studies in crystalline^-^-^-' and thin film,'"" samples 
confirmed the existence of the mid-infrared (MIR) excitation 
for a large number of manganites. In this context the concept 
of small (SP) and large polarons (LP) as essential ingredients 
to understand the electronic transport properties and the MIR 
excitation (labeled polaronic peak in the following) became 
subject of considerable theoretical effortsAi2ii^*ii In contrast 
to an optical transition between electronic bands, a character- 
istic feature of a polaronic absorption process is the significant 
frequency shift with temperature (see, for example, Ref. ll5l) . 

The shape of a polaronic resonance in the optical conduc- 



tivity is predicted by theory whereby the details strongly de- 
pend on the nature of the polaronic excitations: large po- 
larons exhibit a threshold behavior and a coupling dependent 
structure '^''^ whereas small polarons display an asymmetric 
Gaussian peak;'^''' these latter polarons may even pair to form 
bipolarons with their own characteristics.^'*'^" The assignment 
"large" or "small" polaron usually refers to the spacial exten- 
sion of the lattice distortion induced by a single charge carrier. 
True to the literature (for example, see the review by Emin 
on optical properties of polarons-') one would denote a lattice 
distortion generated by an electronic charge as a small polaron 
if it collapses to a single site — the smallest state compatible 
with the atomicity of a solid. Correspondingly, all polaronic 
states with extension larger than a lattice spacing are labelled 
large polarons. A large polaron does not necessarily extend 
over a large number of sites. However, since we have no di- 
rect means of determining the size of the polaronic cloud from 
optical spectroscopy, we instead refer to the related shape of 
the absorption spectra: threshold behavior for large polarons 
versus asymmetric Gaussian peak for small polarons.^' 

For the metallic manganites, a clear identification of the po- 
laron type and the respective model is disputable not only be- 
cause other electronic correlations may correct or even replace 
the standard polaronic physics but also because the statistics 
of the rather dense polaron system may strongly influence the 
shape of the resonance. Actually, small polarons which are 
trapped in a locally polarized lattice are less affected by fi- 
nite doping. The MIR absorption may still be visualized in 
terms of transitions between the adiabatic levels of neighbor- 
ing sites. The presence of other localized polarons at more 
distant sites will renormalize the polaronic energy scale (the 
binding energy) and it will raise the high-frequency wing 
of the resonance due to additional multi-polaron processes. 
However, a finite density should affect less the low-frequency 
wing, that is, absorption at energy less than the binding en- 
ergy. In this work we will provide further optical conductivity 
and dc -resistivity data which confirm this small polaron see- 



nario for thin films of La2/3Cai/3Mn03 (LCMO). 

By contrast, a large-polaron resonance should be influenced 
profoundly by a finite density of polarons due to their spacial 
extension and possible delocalization. A perturbational ap- 
proach for small electron-phonon coupling at finite charge car- 
rier density was presented recently by Tempere and Devreese 
(TD)iS In fact, if we compare the one-polaron absorption of 
Gurevich, Lang and Firsov (GLF)"" to the finite density re- 
sult of TD, we notice that the low frequency behavior above 
the threshold is dictated by different power laws. Whereas 
the one-phonon absorption is characterized by a square root 
behavior, the fermionic statistics of many polarons leads to a 
linear density of states of particle-hole excitations at the Fermi 
edge which is in fact reflected in the threshold behavior of the 
absorption. In a recent paper— we proposed to identify the 
MIR resonance in thin films of La2/3Sri/3Mn03 (LSMO) by 
large polarons of this type. 

Here we extend the preliminary discussion of our previ- 
ous work in order to elucidate the distinct polaronic behav- 
ior in the metallic phase of manganites with strong (LCMO) 
and intermediate-to- weak electron-phonon coupling (LSMO). 
The purpose of our discussion is also, to a certain extent, a 
critical review of those polaronic models which have been 
proposed to identify the resonance in previous measurements 
with bulk samples. We observe that hitherto the proposed the- 
oretical models are not capable of capturing all the aspects of a 
polaronic metal like LSMO. Notable details of the MIR spec- 
tra are in disagreement with the available theoretical models. 
A more comprehensive theory which offers an explanation 
particularly for the temperature-dependent shift of the pola- 
ronic resonance and for the magnitude of the spectral weight 
on the high-energy side would probably clarify the nature of 
the polaronic metal. 

Using thin films in contrast to single crystals, we achieve a 
lower conductivity due to grain boundaries and strain^ Con- 
sequently screening is effectively reduced — to a level where 
phonons and polarons are well observable in the optical con- 
ductivity, even in the metallic phase. This allows the detailed 
analysis of phonons^^'^'' and polarons^^ from the respective 
optical data. Besides, the optical properties of single-crystal 
samples are more sensitive to surface preparation which, in 
particular, influences the MIR-regioni2i 

The role of the lattice mismatch on the structural and trans- 
port behavior in dependence of the film thickness has been 
systematically investigated in numerous articles. In particu- 
lar it depends on the substrate/film combination and on the 
growth conditions. In general thicker films show a behavior 
very close to that of bulk— whereas substrate induced effects 
become more important in thinner filmsi^ 

The structural differences between La2/3Cai/3Mn03, 
which is orthorhombically distorted, and La2/3Sri/3Mn03 
with its rhombohedral symmetry represent promising condi- 
tions to distinguish between competing polaronic models. We 
discuss those models which adequately reproduce the low en- 
ergy side of the resonance, as these excitations characterize 
the nature of the observed polarons. 



II. EXPERIMENTAL ANALYSIS 
A. Details and Characterization of the Films 

The thin films were grown by a standard pulsed laser de- 
position technique^'^ on plane-parallel 10 x 10 mm^ single 
crystafline substrates: LSMO on (LaA103)o.3(Sr2AlTa05)o.7 
(LSAT) and LCMO on NdGa03. Thin films with thickness 
between 200 nm and 400 nm were investigated. The nomi- 
nale lattice mismatch is <0.1% for both. 

X-ray diffraction (XRD) measurements were performed us- 
ing a Bruker D8 Discover four circle x-ray diffractometer. In 
Fig. [2 we show x-ray diffraction patterns (9 — 20) obtained 
for LCMO on NdGa03 (upper panel) and LSMO on LSAT 
(lower panel). For LCM0/NdGa03 the XRD displays only 
(lOl) diffraction peaks for Pnma and for LSMO/LSAT only 
(/OO) diffraction peaks for Pm3m symmetry. For both com- 
pounds there exist no indications of any impurity phase. This 
points to the highly oriented crystal growth of the films. The 
insets of Fig. [2 exhibit the rocking curves of the (202) peak 
for LCMO and the (200) peak for LSMO. Both films show 
a very narrow full width at half maximum (FWHM) of the 
rocking curve of 0.04° which give evidence for a high quality 
crystalline structure of the films. 

Reflectivity measurements of the film-substrate system 
were carried out in the range of 50 to 4 x 10'* cm~* using the 
Fourier transform spectrometers Bruker IFS 1 13v and Bruker 
IFS 66v/S. Below room temperature the measurements were 
performed with a He-cryostat where the sample is placed in 
a static He exchange gas. The spectra at higher temperatures 
were taken in a home made oven in which the sample was ex- 
posed to a continuous flow of heated nitrogen gas. In addition, 
a Mach-Zehnder interferometer was used for ac measurements 
between 10 and 30 cm^* which allowed to obtain both the in- 
tensity and the phase shift of the wave transmitted through 
the films on the substrate. Such transmission measurements 
of thin films have proven to be a powerful method for study- 
ing the electrodynamic properties of highly conductive thin 
films at low frequencies.^' The combined data sets improve 
the Kramers-Kronig (KK) analysis of the film-substrate sys- 
tem in the complete frequency range. For the analysis, the re- 
flectivity R{u!) was extrapolated by a constant for frequencies 
below 10 cm^'^. For the high frequency region, the reflectiv- 
ity between 4 x 10"* cm^'^ and 10^ cm^'^ was extrapolated by 



a uj 



-1.5 



law. Above, a u;^ dependence was assumed. The 
same procedure was utilized in separate experiments for the 
respective substrate. Applying the Fresnel optical coefficients 
for a two-layer system, the optical conductivity of the film was 
evaluated. A more detailed description of the data analysis is 
given in Ref . ^, 

For room temperature the combined data sets of IR and ac 
reflectivity R are shown for the film-substrate system and for 
the bare substrate in the upper panel of Fig.|2]while the calcu- 
lated conductivity a is displayed in the lower panel. Most 
phonon resonances of the substrate are located in the min- 
ima of the film-substrate spectra. In the metallic phase the 
screening shields the lattice modes of the substrate, therefore 
only resonances of the thin film are visible. The temperature 
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FIG. 1 : (Color online) Room temperature — 20 XRD pattern for 
LCMO on NdGaOa (upper panel) and for LSMO on LSAT (lower 
panel). The insets show the rocking curve of the (202) peak and the 
(200) peak, respectively. 



dependence of the phonon resonances provides for an exact 
determination of the film properties independently of the in- 
fluence of the substrate. Moreover, the resonances known 
from neutron scattering measurements''-' are in good agree- 
ment with our results. 

The electrical dc -resistivity of the films was measured with 
a standard four-point technique. ^^ We estimate the defect level 
in relation to the single crystal. Due to stronger strain effects 
and a higher density of grain boundaries, the defect level in 
LSMO #2 is about 12 times higher and in LSMO #1 about 24 
times than in the single crystal. For LCMO a single crystal 
was not available but the residual resistant is almost the same 
as for LSMO #2. 

The magnetization was obtained with a commercial Quan- 
tum Design SQUID magnetometer The measurements were 
carried out for the substrate and for the deposited film. The 
data were reported in a previous publication,^* and a single 
crystal measurement is found in Ref. 25 for comparison. To 
characterize the ferromagnetic to paramagnetic (FM-PM) and 
the metal-insulator (MI) transition of the samples, we define 
Tc as the magnetic Curie temperature and Tmi is determined 
from the maximum of the dc resistivity curves. For LSMO the 
onset of spontaneous magnetization appears at 345 K, well be- 
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FIG. 2: (Color online) Upper panel: Reflectivity of the film-substrate 
system (LSMO on LSAT) and of the substrate itself (LSAT) at 295 K. 
Lower panel: The optical conductivity a of the film, evaluated from 
the reflectivity data. 
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B. Experimental Results 

The reflectivity spectra for frequencies below the optical 
transitions are displayed in Fig.|3]for LSMO (upper panel) — 
and for LCMO (lower panel). The solid lines represent the 
spectra for temperatures in the metallic FM pase while the 
broken lines refer to those in the PM phase. The series of 
sharp peaks below ~ 700 cm~^ were identified as phonon 
modes. A detailed analysis of the phonon excitations in these 
films has been presented in a previous articled The identi- 
fication of the infrared-active phonons attests the quality of 
the thin films. In the following we focus on the excitation 
above ~ 700 cm^^, which has been assigned to polarons 

,..f..j.^ 2.3.6.7.8.9.10.11.23 



For increasing temperatures the polaronic resonance in the 
optical conductivity continuously shifts to higher frequencies 
and, moreover, a loss of spectral weight is observed — as has 
been discussed in Ref. 23. Both quantities, frequency and am- 
plitude of the resonance, will be listed in Tabs. HI and IIII for 
LCMO and LSMO, respectively, linked to specific theoreti- 
cal models. There is a remarkable difference in the shape of 
the polaron peak for LSMO and LCMO. In LSMO the shape 
is asymmetric and displays a steep slope below the peak posi- 
tion, whereas for LCMO the broader hump is more symmetric 
about its maximum and a low frequency threshold is entirely 
missing. 

When we scale a{i') by a typical energy i^o of the polaronic 
resonance we find a universal low-energy slope for the LSMO 
spectra (upper panel of Fig.|4} whereas this scaling behavior 
is absent in the LCMO spectra (Fig.|3 lower panel). For this 
presentation, i^o is the threshold on the low-energy side of the 
MIR-resonance of the LSMO spectra and the frequency at the 
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FIG. 3: (Color online) Reflectivity spectra for LSMO on LSAT (up- 
per panel) and LCMO on NdGaOs (lower panel) over a temperature 
range from the metallic to the insulating phase. The solid lines rep- 
resent spectra in the FM phase, the broken lines in the PM phase. 



FIG. 4: (Color online) Rescaled optical conductivity of LSMO (up- 
per panel) and LCMO (lower panel): the optical conductivity is 
scaled by the maximum of the MIR resonance. A constant back- 
ground of (Tbg =370 n~^cm~^ for Act is subtracted for LSMO (for 
details, see main text). The frequency is scaled by the threshold value 
f of the onset of the polaronic peak for LSMO, while for LCMO by 
the peak maximum. Clearly, the scaling behavior of the low fre- 
quency slope is missing for the LCMO spectra. 



III. FURTHER EVIDENCE FOR SMALL AND LARGE 
POLARON SCENARIOS 



maximum of the MIR resonance of the LCMO spectra, re- 
spectively. To retrieve the threshold frequency fo for LSMO, 
one has to take into account that a background adds to the 
polaronic absorption close to the threshold. We assume that 
the 295 K data approximately display the background for the 
considered frequency window: the phonon resonances are suf- 
ficiently far below j/q so that they should interfere very little 
with the polaronic resonance, and the frequency dependence 
of the intensity of excitations responsible for the background 
is negligible for the considered frequency window in which 
we address the scaling behavior We then extrapolate the low 
frequency slope of the curves at all other temperatures to this 
"zero line" in order to identify uq. Strictly speaking, i^o de- 
pends slightly on the extrapolation procedure. This qualita- 
tive distinction of the LSMO and LCMO spectra has to be 
addressed in a thorough theoretical modelling. 



A. Holstein SP hopping 

The MIR optical conductivity displays distinct polaronic 
excitations for LSMO and LCMO. It is expected that this dis- 
tinction also manifests itself in a disparate dc-resistivity for 
the two types of manganite films. 

In Fig. Is] we show the high temperature dc-resistivity of 
LSMO (upper panel) and LCMO (lower panel) in a lii{p/T) 
vs. 1/T scale. In the lower plot a linear dependence for 
the latter compound above Tmi is clearly visible. Slightly 
above Tmi, which coincides with the magnetic ordering tem- 
perature, deviations from the thermally activated behavior oc- 
cur It is widely accepted that the electron-phonon coupling is 
sufficiently large to localize the carriers in the paramagnetic 
phase — presumably in combination with other localization 
mechanisms. Over a temperature range from 255 — 560 K 



the electrical dc-resistivity is well described by the adiabatic 
SP model^^ (lower panel) with a dc activation energy Ea = 
672 cm~^, which is in agreement with previously reported 
results. ^"•^"^'^''-^^ In general the binding energy Ei, of the SP in 
the optical data (see Tab. EJi is related to the activation energy 
of the transport data by Eb = 2Ea- The maximum in the op- 
tical conductivity in the paramagnetic phase (T — 250 K) is 
located around 6040 cm^^, which is rather 8Ea than the the- 
oretically expected 2i5b = AEa- The discrepancy between the 
transport value and the optical value is not surprising since it 
was also observed in other polaronic systems'*** (including the 
classical small polaron material Ti02; the difference has been 
discussed by Reik^'' and Austin^ and for the manganites by 
Hohenadler and Edwardsii). 

Neutron scattering measurements for single crystal sample 
were interpreted in terms of "sin gle polarons" and "correlated 
polarons" (see, e.g. Ref. I42l43l) — the existence of which is 
supposed to be consistent with a phase separation scenariojii 
Alvarez et al^^ predict uncorrected magnetic clusters in a 
temperature range Tc < T < T*. Below Tc the clusters 
align their moments. Above T* single polarons are expected 
to dominate the transports" In fact, for LCMO we find a kink 
in the high temperature dc-resistivity data (see In p{T) in the 
inset of Fig. |5}. We tentatively identify this kink at 560 K 
with T*. This observation is qualitatively consistent with film 
and single crystal results reported in Ref.|2Jbut we question 
their interpretation according to which the kink is a signature 
of a structural transition. Rather it appears that the different 
slopes depend on various scattering processes of correlated 
and single polarons, and the temperature T* characterizes the 
crossover. 
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Previously, it has been reported that correlated polarons 
have been found only in the orthorhombic phase (LCMO).^— 
This observation may be related to the distinct high tempera- 
ture behavior in LCMO and LSMO. For LSMO we identified 
LP-type excitations from the optical conductivity data. Fig- 
ure|5](upper panel) shows the high temperature dc-resistivity 
from 330 — 750 K. Clearly, the adiabatic SP model (solid line) 
provides a poor fit to the experimental data (compare the dis- 
cussion in Ref. 1471). 

The temperature dependence of the JT distortion for 
Lao.ySro.sMnOa has been investigated more recently by Man- 
nella et alJ^i^ with x-ray absorption spectroscopy. They pro- 
vide direct experimental evidence for the presence of polarons 
in the paramagnetic state and analyze the dopant dependence 
of the JT distortion. Specifically, the variance of the Mn-O 
bond length distribution in the Sr doped sample was shown to 
be only half that in the Ca doped sample. 



B. Charge Carrier Mobility 

Further support for the SP concept in LCMO and the LP 
ansatz in LSMO comes from the literature on the polaron mo- 
bility. Figure |6l compares the temperature dependence of the 
dc- and ac-conductivities (right axis) with the mobility (left 
axis) obtained by Hall-effect measurements by Jakob et al^ 
The carrier mobility decreases with increasing temperature in 
the FM phase. Above Tmi the mobility increases slightly. All 
three measurements show the same overall temperature be- 
havior. The deviation close to Tmi may be due to different 
coupling mechanisms of the local transport properties to crit- 
ical magnetic fluctuations. 

The decrease of the mobility in the FM phase corresponds 
to the increase of the effective electron-phonon coupling 
(polaron energy over bandwidth). Moreover, the spectral 
weight of the polaronic excitation is reduced on increasing 
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FIG. 5: (Color online) High temperature dc-resistivity of LSMO (up- 
per panel) and LCMO (lower panel). The solid lines are the fit of the 
adiabatic SP model. For LCMO the curve is linear and corresponds 
to polaron conduction with an activation energy Ea = 672 cm^^. 
The inset displays a kink in the dc-resistivity of LCMO at 560 K 
which we identify with T* (see main text). 



FIG. 6: (Color online) For LCMO the temperature behavior of dc and 
ac electrical conductivity a (right axis) corresponds to the tempera- 
ture dependence of the mobility (left axis) obtained from Hall-effect 



measurements." 
ical origin. 



The similarity of the data indicates a common phys- 



temperature. As a consequence, with decreasing lattice dis- 
tortions, the tendency to charge locaHzation is reduced. 

For LSMO only few Hall effect studies have been re- 
ported, none of them providing values of the polaron mobility 
/i/f ."^^-^'-^^ Therefore, only a rough estimate as Rh / Rxx ~ 
400 mm^/Vs at 4.2 K is possible^^ which is four times larger 
than in LCMO (Fig. |6j. The ratio of the mobilities corre- 
sponds well to the ratio of the effective masses of specific- 
heat measurements,^^ and reveals the expected tendency. We 
propose that optical and Hall transport consistently support 
the applicability of a LP-type approach. This interpretation is 
also in agreement with neutron scattering results where the 
spacial extension of the polaron covers several lattice con- 
stants at low temperatures: Louca and Egami^'* investigated 
deviations of the local structure from the average crystal struc- 
ture in Lai_a;Sra;Mn03 by pulsed neutron scattering with pair 
distribution function analysis. X-ray absorption fine structure 
measurements^^'^^ differ for the metallic phase quantitatively, 
but all groups identify a temperature-dependent lattice distor- 
tion. 



IV. DISCUSSION OF THE MIR SPECTRA 

Early on the double exchange (DE)piS in combination with 
a strong Hund coupling, was identified as a key element to 
explain the interrelation between metallicity and ferromag- 
netism in the doped manganites. It has only recently been 
understood that further basic mechanisms have to be an inte- 
gral part of a comprehensive model that can explain magne- 
totransport properties and the transition to the paramagnetic 
phase. Millis, Littlewood and Shraiman^^ elaborated that the 
DE mechanism is not sufficient to reproduce the measured re- 
sistivity for temperatures of the order of Tq. They suggested 
that, in the presence of a strong electron-phonon coupling, 
polaronic effects have to be included. Dagotto, Hotta and 
Morec^^ promoted that intrinsic inhomogeneities are gener- 
ated through the cooperative effect of phase separation ten- 
dencies, long-range Coulomb interaction and disorder, and 
these inhomogeneities are supposed to dominate the transport 
behavior 

In the following discussion on the nature of the MIR reso- 
nance, we will consider only the polaronic models with no re- 
finement of electron-electron interactions or disorder physics. 
This disregard of some of the essential ingredients of a com- 
prehensive microscopic model for manganites is in general 
not justified. However, a solution to such a "comprehensive 
model" is at present not available, certainly not for an authori- 
tative evaluation of the MIR optical conductivity. On the other 
hand, some of the electronic correlations missing in the po- 
laron models would enter only as effective parameters, such 
as the temperature-dependent kinetic energy, or the projection 
on one spin direction. Consequently, these correlations do not 
necessarily have to be disentangled from a more microscopic 
modelling. Finally, the purpose of our discussion is also a 
critical review of the polaronic models which have been pre- 
viously related to the resonance. We reemphasize that in the 
considered thin films, the resonance is more pronounced and 
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FIG. 7: (Color online) The optical conductivity spectra of LSMO 
(upper panel) and LCMO (lower panel) at 100 K (data points) on 
a logarithmic frequency scale. For a preliminary assignment to 
large polarons (LP, solid line) and small polarons (SP, dashed line), 
Eqs. jsj and Q have been used for the respective fits. 



the shape is well resolved so that a refined analysis is now 
attainable with our data. 

A small polaron in the regime of sufficiently strong 
electron-lattice coupling is excited from a site-localized state 
to another well-localized state at an adjacent site.^^ The max- 
imum of the absorption, that is the most probable process, 
is the Franck-Condon like transition with frequency LOm four 
times the activation energy for the intersite transition. ■''^•^" The 
MIR conductivity then has the form of an asymmetric Gaus- 
sian peak centered around cj„i. Whereas the optical conduc- 
tivity of LCMO approaches this shape and therefore is iden- 
tified as a small polaron, the optical conductivity of LSMO 
does not (compare Fig.Q. 

A large polaron is observed usually in the weak coupling 
regime. Then, to first order in the electron-lattice coupling, 
the optical absorption is the excitation of a phonon with a si- 
multaneous excitation of the charge carrier GLF'^. Conse- 
quently, one finds a threshold behavior at the (optical) phonon 
frequency with a sharp increase in the absorption rate due 
to the available scattering states and a slow decrease at the 
high-energy side controlled by a l/o^'^-factor, most easily 
identified in the force-force correlation function form of the 



conductivity.^" For intermediate to strong electron-lattice cou- 
pling, if the polaron is in a self-trapped state, the threshold is 
determined by the polaron binding energy, with the charge 
carrier scattering into an extended state (Emin^' assumes a 
free-carrier state). The initial upturn of the absorption spec- 
trum is again given by the phase space of the available scat- 
tering states. The shape of the observed optical conductiv- 
ity appears to be qualitatively similar to that in the GLF- and 
Emin-approaches, and we tentatively identify the polaronic 
resonance for LSMO as the excitation of a large polaron (see 
Fig. 0. The validity of such models to the manganite films 
will be investigated in the following sections. 

The asymmetric shape of the MIR resonance with the steep 
increase on the low-energy side is quite spectacular for the 
investigated LSMO films. As argued before, the upturn of 
the optical conductivity above the threshold is controlled by 
phase space considerations. For example, if the electronic 
states created by the absorption process are particle-hole ex- 
citations close to a Fermi edge, the number of such states in- 
creases linearly with the energetic distance (frequency) to the 
edge. Whatever the details, one may expect a scaling be- 
havior of LSMO, at least for the initial increase above the 
threshold. This scaling would signify that the low energy scat- 
tering processes are of the same origin, independent of tem- 
perature, related band narrowing and other energy scales. In 
fact, a universal low-frequency slope has been verified for the 
LSMO specti-a (Fig.0] upper panel) but it is absent for LCMO 
which neither displays a threshold on the low-energy side of 
the MIR-resonance ( Fig.|2] lower panel). 



A. Small polarons in LCMO 

Various experimental data sets on both, ferromagnetic and 
paramagnetic manganites, have been interpreted in terms of 
the Holstein small-polaron theory^^ (see our previous section, 
and also the critical discussion on this issue in Ref. 611). An 
extension of the standard SP optical conductivity formulaiSiS 
has been proposed for the evaluation of the low-temperature 
optical spectroscopyi^2iS 
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A = 2V2^6^vib 

Here Ei, is the SP binding energy, ct(0, T) is the dc con- 
ductivity, A signifies a broadening, and i?vib is the char- 
acteristic vibrational energy, which is i?vib = fc^T in the 
high-temperature regime and E'vib = ?iijJph/2 at low tem- 
peratures {ksT < huiph). Eq. constitutes a broad reso- 
nance with a maximum at hujm sa 2Eb. For high tempera- 
tures, the dc-conductivity results from the thermal activation 
over the potential barrier {Eh/2) between adjacent sites and 
the optical conductivity peaks at the frequency {'2,Ei,) of the 
Franck-Condon-like transition for neighboring potential en- 
ergy curves, which are approximately parabolic in the con- 
figurational coordinates. For low temperatures, the processes 



are similar, however the phonons which participate in the ac- 
tivated hopping are not thermal. 

The dashed lines in Fig. present the optical conductiv- 
ity of a SP, Eq. Q. The shape for low frequencies, includ- 
ing the maximum, is well reproduced for LCMO. The SP 
model nicely describes the experimental results for 100 K us- 
ing 2Eb = 2820 cm^^ and a characteristic phonon energy 
£'vib = 270 cm^^. The disagreement between theory and ex- 
periment for higher frequencies is not unexpected (it has also 
been observed for the classical small polaron material Ti02, 
and several mechanisms were proposed to account for the 
difference^ see also the SP evaluation for the 2d nickelates 
and cuprates-^). As seen from Tab. U and Tab. |lll| with ris- 
ing temperature the polaron formation energy increases. For 
LCMO the binding energy i2Eb) rises from 2210 cm^^ (6 K) 
to 6040 cm^ ^ (250 K). Above the metal-to-insulator transition 
at Tmi, the polaron formation energy remains constant. 

Whereas one may expect that the low-energy SP absorp- 
tion is well explained by the processes described above — 
since they are local — , one has to consider the consequences 
of a high polaron concentration and Fermi energy for the 
high-energy SP absorption in order to understand the slow 
decay observed in experiments. Especially for temperatures 
of the order of Tc polaron-polaron interaction appears to be 
strong,"*^"*^ however a relation to the high-energy slope of the 
polaronic MIR resonance is not evident. 

An explanation of the high energy decay and of the tem- 
perature dependence of the polaron binding energy is beyond 
the scope of this paper A transfer of spectral weight to- 
wards higher energies was discussed within a dynamical mean 
field theory of the Holstein model by Fratini et al.,^^ how- 
ever, their treatment only applies to a dilute nondegenerate 
polaron gas. The shape and position of the polaronic MIR 
resonance for a model with Fermi edge and a finite concentra- 
tion of polarons has been investigated by Millis, Mueller and 
Shraiman (MMS)."* They study the DE model with a JT form 
for the electron-phonon coupling. The local phonons are taken 
as classical harmonic oscillators, and the electronic correla- 
tions are evaluated within the DMFT scheme. For half filling 
(rie = 1, with the Cg orbitals singly occupied) and intermedi- 
ate coupling, MMS compare their optical conductivity favor- 
ably with the data of Okimoto et al?-^^ for Lao.sysSr lysMnOs 
— and for somewhat stronger coupling with the data of Ka- 
plan et al? for Ndo.ySr sMnOa. Also our optical data agree 
qualitatively with the latter results of MMS though the high 
frequency side again decays too fast in the model evaluation. 
The broad peak is due to transitions between the two JT split 
levels with a strongly fluctuating phonon coordinate. The ac- 



TABLE I: Fit parameters of the optical conductivity for LCMO at 
different temperature using the SP formula, Eq.Q 
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tual results for the evaluation of doped systems are less favor- 
able: a second peak below is formed which is due to transi- 
tions to unoccupied sites with no JT splitting. It appears to 
contradict the observations. However, MMS argue that fur- 
ther electronic correlations which have not been included can 
remedy this shortcoming. 

Some of the quantum fluctuations missing in MMS are in- 
cluded in the "many-body CPA" treatment of the Holstein- 
DE model^^ for which Hohenadler and Edwardsii compare 
the MIR optical conductivity with the measurements^ for 
Ndo.ySr.sMnOs. The shift of spectral weight to lower fre- 
quencies with decreasing temperature is still insufficient. Per- 
roni et al.^^ propose that their variational Lang-Firsov treat- 
ment recovers the low temperature behavior,^'''^** including a 
Drude peak. 



B. Bipolaron theory for the optical conductivity 

The crossover from the paramagnetic to the ferromagnetic 
phase in the optical absorption has been studied by Alexan- 
drov and Bratkovsky'"* in terms of a bipolaron model in 
which the optical intraband conductivity should be provided 
by two different polaronic contributions originating from SP 
and bipolarons with different binding energies. Accordingly, 
in the PM phase the optical conductivity would be dominated 
by small bipolaronic charge carriers, which break up at the 
magnetic transition, leaving only SP which then dominate the 
MIR excitations in the FM regime. As a result, a sudden trans- 
fer of spectral weight at Tc is predicted along with a sharp 
drop in the dc -resistivity due to the charge carrier collapse. 
As such a transfer is not confirmed by our data, we rather dis- 
card such a scenario for LSMO and LCMO, even though the 
oxygen-isotope effect on the thermoelectrical power in sev- 
eral ferromagnetic manganites seems to provide experimental 
clues for bipolaron formationj2ii 



C. Large polarons in LSMO 

We will now inquire if models, which have been recently 
proposed in this context, indeed reproduce the observed scal- 
ing and threshold behavior for LSMO. For this purpose, we 
organize this Section lTV Cl as follows: 

In a first approach we will use a well-controlled weak cou- 
pling scheme (see Sec. IIVC 1> which follows from a RPA- 
evaluation of the Hubbard model to which a Frohlich type 
electron-phonon coupling has been added. As we only con- 
sider the ferromagnetic metallic state we restrict the evalua- 
tion to the spin polarized sector (spinless fermions). The local 
interaction U signifies the on-site interaction between charge 
carriers in the two e^-orbitals. The RPA-results well approxi- 
mate the low-frequency slope of the polaronic resonance and 
reproduce the observed scaling. In fact, the model replicates 
the most prominent and important characteristics of the po- 
laronic spectra in LSMO. However it fails to treat the high 
energy features adequately. 



In a second approach we investigate intermediate and 
strong coupling schemes for large polarons (see the end of 
Sec. lIVCTI and the Section irVC2l on self-trapped polarons). 
The latter scheme was proposed recently,^ in order to identify 
the MIR absorption processes in bulk manganites as excita- 
tions from a localized self-trapped state to free-carrier contin- 
uum states (the authors of Ref. 7 refer to the model presented 
by Emin in Ref. 21). The shape of the polaronic resonance 
in LSMO is well reproduced for the low energy side and the 
model may account for a shift of the threshold energy with in- 
creased binding (at higher temperatures). However the theo- 
retical (strong coupling) estimate of the polaron radius differs 
significantly from the radius which fits the experimental data. 
Moreover, the approach does not apply when a Fermi surface 
of the (free) charge carriers is formed. 

Finally, in Sec. lIVC"3] we discuss the consequences when 
charge carriers are excited from a potential well to a Fermi 
surface. The collective behavior associated with the response 
of the Fermi sea to the transition determines the shape of the 
resonance in the adiabatic limit (in which the potential well 
decays slowly with respect to the fast electronic transition). A 
full evaluation of this fascinating scenario is not yet available. 



1. Weak coupling scenario 

The shape of the MIR resonance in the LSMO films is rem- 
iniscent of the weak couplin g res ult by GLF for Frohlich po- 
larons (see, for example, Ref.l60land Eq. @ below). A closer 
inspection reveals that the width of the GLF absorption is too 
small and the low-frequency slope is too steep to fit the LSMO 
data (see Fig.|8}. This failure is not unexpected since we have 
10^^-10^^ charge carriers per cm"^, as identified from Hall 
measurements.^' An extension of the GLF weak-coupling po- 
laron theory to a finite concentration has been carried out by 
Tempere and Devreese.^^ In their work on optical absorption 
of an interacting many-polaron gas, they found a remarkable 
agreement between their theoretical prediction and the exper- 
imental curve for Nd2Cu04_j,^' in the optical conductivity 
around 1000 cm^'. We now transfer their approach to the 
manganites. 

The real part of the optical conductivity <j{uj) is expressed 
through a momentum integral over the dynamical structure 
factor S{q, lo) which is in three space dimensions i^2iffi 
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Here, Up is the polaron density, m is the effective mass of the 
charge carriers, a is the electron-lattice coupling and hujo — 
hvQ is the threshold energy for polaronic absorption (in weak 
coupling, it is the energy of the respective optical phonon). 
The dynamical structure factor is expressed in terms of the 
dielectric constant e{q, lo) 
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where n is the density of charge carriers which should be iden- 
tified with the polaron density rip for the considered finite con- 
centration. TD evaluated Eq. (|2j to linear order in a which im- 
plies that S{q,uj)is taken to zeroth order in the electron-lattice 
coupling. The dielectric constant is then approximated by the 
Lindhard (or RPA) expression for the homogeneous electron 
gas. In this linear approximation (with respect to a), TD ar- 
gue that, for the low-density limit, the GLF result is recov- 
ered. Indeed, for fermions with one flavor, the structure factor 
approaches'^ S{q,uj)kf^o ^ 2nS{oj — Eq/Ti). The optical 
conductivity assumes the limiting form of GLF"' (cf. Figlsjl: 
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For the manganites in the investigated ferromagnetic, metal- 
lic regime with a carrier density of 10^^-10^^ cm^'^, we take 
the following standard simplifications for the low- temperature 
evaluation: (i) only one spin direction prevails (spinless 
fermions), due to DE and strong Hund coupling, (ii) there 
are two degenerate orbital states (the two e^-levels); accord- 
ingly we deal with two flavors for the fermionic particles, (iii) 
the interaction is represented by the local Hubbard U for two 
fermions on the same site in two tg orbital states, and not by 
the long-range Coulomb interaction Vq = Ane'^/q^, (iv) the 
fermions acquire a quadratic dispersion. The last assumption 
may be refined to the appropriate two-band dispersion of the 
Bg-levels. However, for the discussion of the optical conduc- 
tivity it is of minor importance, especially when we primarily 
discuss the low-frequency excitations with energy much less 
than the Fermi energy. Furthermore, a dynamical JT effect is 
not implemented: it is entirely missing in this weak coupling 
approach of TD. 

The longitudinal dielectric response in the Hubbard model 
(with quadratic dispersion) is in RPA: 
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The related expression for the Coulomb gas'^ is usually not 
presented in a closed form. Here care has been taken to ensure 
that this closed expression for the dielectric response is valid 
for all complex values of y. 

The optical conductivity is (cf. the similar expression in 
Ref.|23for the Coulomb gas model) 
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FIG. 8: (Color online) Comparison of the low temperature (100 K) 
optical conductivity (in the MIR range of the polaronic excitations) 
to the results from various model calculations: the continuous line 
refers to the weak coupling approach of Tempere and Devreese mod- 
ified for an on-site Hubbard interaction, Eq. J6}; the dashed line is 
the result of the phenomenological approach for self-trapped large 
polarons by Emin, modified to Gaussian localized states, Eq. J9}; the 
dotted curve is the weak coupling single-polaron result by Gurevich, 
Lang and Firsov, Eq. J4j. The appropriate parameters for each of 
these model evaluations are given in the main text. 



whereby 
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3 m ojQ 
is the prefactor which already appears in the GLF-limit, and 
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is a dimensionless prefactor which accounts for the finite 
charge carrier density. The integration variable x in Eq. (|6j is 
the normalized wavevector a; — q/kf. The dielectric response 
in the Hubbard model, Eq. (|5}, and the optical conductivity, 
Eq. (|6j, are in a simple way related to the respective expres- 
sions for the Coulomb gas if we substitute in Eqs. (|5} and (|6j 
the interaction parameter by 
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where ^ is the dimensionless interaction strength of the 
Coulomb gas, ^ = e^kf /{-KEf). 

Finally, we adjust the parameters tentatively in the follow- 
ing way: a) the density of charge carriers is 6 x 10^^ cm^^, 
which is the stoichiometric number of doped holes and 
which is in the range of what has been estimated from Hall 
measurements,^' b) we take m three times the bare elec- 
tronic mass (m — 3me) which is consistent with specific heat 
measurementspi c) since Ef = 3270 cm^^ for the density of 
(a) and for the effective mass of (b), and since we read from 
the 100 K data i^o = 617 cm^^ we fix hujo/Ef = 0.19, that 
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is, we are in the adiabatic regime, and d) U/Ef — 0.2 in order 
to start, for a first estimate, with small electronic interactions 
for which RPA is a reasonable approach. 

For this parameter set we find from Eq. ^ a shape of the 
polaronic resonance (see Fig.|S) which reproduces well the ex- 
perimental data for frequencies up to the maximum. The high 
energy slope, however, decreases more rapidly than in the 
measured optical conductivity. The maximal optical conduc- 
tivity is, in our calculation, Umax = a x 1.84 x lO"^ fl^^cm^^ 
which should fit the data with a ~ 1.3 — whereby the shaded 
background in Fig.|8lis to be subtracted. 

If we had taken the Coulomb gas (as TD did for 
Nd2Cu04-y) with the same density, effective mass and the 
same ratio Tiuoo/Ef as above for the Hubbard model, we 
would find nearly the same line shape (not distinguishable 
in the plot). The dimensionless interaction strength then is 
^ = 6.5 and the plasma frequency is at fioj plasma/ Ef — 4.15. 
Since the plasmon is consequently far outside the frequency 
window of our investigation it is not surprising to have similar 
results for the Coulomb-gas and Hubbard-model approach for 
these reasonably large densities. However, the maximal opti- 



cal conductivity is a„ 



a X 0.11 X lO"* J7"^cm"^ in the 



Coulomb gas case. This impUes an electron-lattice coupling 
a of the order of 20, which is probably too large^^ by a factor 
of 5 to 6 — apart from the fact that the Coulomb gas model is 
not quite appropriate for these manganites. We conclude that 
the Hubbard model, rather than the Coulomb gas approach, is 
the proper model Hamiltonian provided that also an electron- 
phonon coupling term is added. 

The solution, Eq. (|5Jl, suggests to determine the model pa- 
rameters through an "optimal fit", now including the high fre- 
quency slope. The shape of the optical conductivity in this 
approach only depends on huso/Ef, and the fine tuning fixes 
the value of U/Ef. For small values of hojo/Ef, such as for 
the above value 0.19, the absorption shape is nearly indepen- 
dent of U/Ef, except that the magnitude (amax) scales down 
with increasing U/Ef 7^ A better fit is found with larger val- 
ues of hujo/Ef. The optimal values are Tiuso/Ef = 0.9 and 
U/Ef ~ 4 (see Fig.|9jl. If we again assume that the density 
of charge carriers is 6 x 10^^ cm^'^, we are constrained to 
m/nif. — 14.3. Alternatively, we could fix m/mf. = 3 but 
then we have to assume a strongly reduced polaron density of 
0.6 X 10^^ cm^'^. The drawback of this "optimal fit" however 
is that (Tniax/a = 0.05 x 10'^ fJ^^cm^^, in the first case, and 
0.02 X 10^ ri~^cm~^, in the second case, both of them too 
small to be reasonably close to experiment. We deduce from 
this investigation that the high energy side cannot be satisfac- 
torily captured within this model. The high energy slope with 
its complex cooperative excitations has to be explored with 
more involved theoretical concepts which necessarily depend 
on a larger parameter set. However, this does not invalidate 
the evaluation for the universal low energy side. 

For the 6 K data, the low energy side of the polaronic res- 
onance was fitted successfully in our previous paper. ^''•^'* For 
the sample presented here we find vq = 593 cm~^, which 
impHes that Tiuo/Ef = 0.18. With U/Ef ^-2 we have 
a nearly identical fit to that of Fig. 4 in Ref. 23 which gen- 
erates sufficient spectral weight so that the electron-phonon 
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FIG. 9: (Color online) In the RPA-approach, Eq. ^, the shape of 
the optical conductivity depends on the model parameters hoJo/Ef 
and U/Ef. For the adiabatic weak-coupling case (U/Ef — 0.2, 
hujo/Ef = 0.19) the approximate form, Eq. jS}, is indistinguishable 
from the full RPA-result, Eq. J6|. For details of the comparison be- 
tween the experimental 100 K data and the theoretical results see the 
main text. 



coupling a can be chosen to be of order 1 . This situation is 
qualitatively similar to that with the 100 K data. The same 
applies to higher temperatures. We may assume that temper- 
atures as high as 295 K are still considerably smaller than the 
Fermi energy. Then temperature does not explicitly enter the 
dynamical response but it is supposed to tune the band width, 
and correspondingly the Fermi energy and effective mass, and 
possibly ujQ (see the discussion below on this last presump- 
tion). 

Assessment of the weak coupling approach. From the above 
discussion we deduce that a value of fiuJo/Ef ~ 0.2 and 
U/Ef < 0.5, is a reasonable choice for low temperatures to 
acquire sufficient optical weight in this approach. Thereby, 
the precise values of U/Ef and Hcoo/Ef just control the 
overall magnitude of the polaronic resonance provided that 
hu!o/Ef ^ 1 applies. ^^ Furthermore, for the considered fre- 
quency range above the threshold, we are always sufficiently 
close to the Fermi edge in order to approximate the frequency 
dependence of the imaginary part of e{q,uj — uiq) hy uj — ujo 
which accounts for the number of excited low-energy particle- 
hole pairs. Moreover, since the frequency dependence of the 
absolute value of e(g, w — uiq) is weak, we expect the optical 
conductivity, Eq. (|6jl, to approach 



a{u) 



aC^ 
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where A is a frequency independent amplitude. Indeed, this 
formula matches the frequency dependence of cr(tj) so closely 
for the considered adiabatic regime of Tiujo/Ef that it is indis- 
tinguishable from the curve generated by the full weak cou- 
pling relation, Eq.|6l(Figs.|8land|9] continuous lines). 

For higher temperatures, Ticuo/Ef <C 1 is not necessarily 
true since u!o{T) increases considerably (see Table HH. Also 
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FIG. 10: (Color online) For 275 K: the RPA-approach, Eq. ^6}, is 
evaluated with a reduced Fermi energy (a reduction factor of 5 as 
compared to the 100 K fit which is the continuous line in Fig.|9). 
The adiabatic weak-coupling form, Eq. jsj' still coincides with the 
full RPA-result and the experimental data along the low frequency 
slope. 



the effective mass of the charge carriers increases with tem- 
perature — due to a reduced double exchange — which im- 
plies a reduced Fermi energy. For the 275 K data we present in 
Fig.llOlthe evaluation of Eq. ^ with a Fermi energy five times 
smaller than for the low temperature fit. The dimensionless 
energy parameters now take the values hcoo/Ef ~ 1.5 and 
U/Ef ~ 1. Obviously, the RPA-result approximates the low 
frequency slope of the polaronic resonance at 275 K with the 
same rigor as for 6 K or 100 K. As anticipated, the adiabatic 
weak-coupling result (Eq. (|8j, see dashed line in Fig. llOt now 
does not exactly match the full evaluation but the deviations 
are still minor To gain the spectral weight of the measured 
polaronic resonance, the electron-coupling has to be larger for 
the 275 K evaluation: a ~ 4 — 5. This is to be expected since, 
with increasing mass, the coupling should also be enhanced. ^^ 
Although the weak coupling approach, Eq. (|S}, is not fully 
applicable for this intermediate-coupling regime at high tem- 
peratures, the result is still so convincing that we propose to 
discuss all data sets, irrespective of temperature, within the 
presented approach. 

The curves for all considered temperatures have the shape 
which is approximated by Eq. (jS) for low frequencies, with 
the threshold uiq replaced by ujo{T). This statement agrees 
with the conjectured scaling behavior for the low frequency 
slope, examined at the beginning of the subsection on large 
polarons — see Fig. 0] in which the low-frequency behavior 
was actually extrapolated to the "zero line" with Eq. (|8}. The 
parameters for this scaling form, the amplitude A{T) and the 
threshold frequency t'o(r), are displayed in Table ITTI 

Shortcomings of the weak coupling approach are the unex- 
plained (temperature dependent) threshold frequency ujo{T) 
and the unmatched slope at the high frequency side when we 
use a reasonable parameter set in order to reproduce the spec- 
tral weight. The coupling parameters are, with this restriction. 



in the small coupling regime. The power law decay at the high 
energy side is u;^'^ in Eq. (|S} whereas the observed power law 
is closer to cj ^ ^ for the 6 K and 1 00 K data (and approximately 
w-0.3 foj. the 275 K data). 

The latter shortcoming should not be taken as a serious flaw 
of the presented concept. The phase space for excitations in- 
creases strongly when the excitation energy deviates from the 
proximity of the Fermi energy and, moreover, the absorption 
through such "high energy" excitations is by no means univer- 
sal. We cannot expect that the simple treatment within RPA 
with quadratic dispersion can still capture these high-energy 
excitations correctly. 

The failure of the weak coupling approach to account for 
the (temperature dependent) threshold frequency ujo{T) is 
more disturbing. The threshold in this approach is fixed at the 
respective frequency of the lattice deformation which will not 
shift from about 593 cm^^ (threshold for 6 K) to 1 183 cm~^ 
(threshold for 295 K). Resolving this issue seems to require 
a type of intermediate-coupling model where the threshold is 
renormalized but the shape of the absorption spectrum is sim- 
ilar to the shape in the weak coupling approach. A possible 
scenario is that of Kartheuser, Evrard and Devreese (KED, 
Ref. Um see also the review by Mitra et al.J'^ section 2.6, 
for further references and for a discussion of the literature on 
this topic). KED propose that the optical absorption for inter- 
mediate coupling consists of a relatively narrow zero-phonon 
peak, assigned to a "relaxed excited state" (RES), and of a 
broad (multi phonon) sideband, with threshold at tJLO above 
the zero-phonon peak. Here, cjlo is the frequency of the LO 
phonon. Consequently the threshold will shift with the fre- 
quency of the RES which is dependent on the electron-lattice 
coupling. This coupling hinges on the electronic band width 
which in turn depends, for the manganites, on temperature. 
This scenario implies that a RES resonance exists at wlo be- 
low the polaronic continuum and that it shifts with tempera- 
ture as the threshold frequency u!o{T) in our data. However, 
such a temperature dependent RES has not been observed as 
an additional sharp resonance for the manganites, neither in 
the optical conductivity nor in Raman scattering. For this rea- 
son we consider the KED scenario inapt to explain the tem- 
perature shift of the threshold in the metallic LSMO films. 

In spite of these shortcomings, it should be appreciated that, 
with respect to the standard polaronic models, this TD-like 
approach implements a Fermi surface. Short-range Coulomb 
interactions are included within a controlled weak coupling 
scheme. Both, Fermi statistics and Coulomb interaction, lead 
to substantial consequences for the shape of the polaronic ab- 
sorption spectrum. 

As we already pointed out, a dynamical JT effect is miss- 



TABLE II: T-dependent parameters of the approximate result for LP, 
Eq. li. 
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ing in this weak coupling approach. Mannella et al. investi- 
gated the temperature-dependent JT distortion (see Fig. 3 of 
Ref. l48h . The temperature dependence of the variance of the 
Mn-O bond length distribution corresponds with the tempera- 
ture dependent shift of the polaronic MIR excitation for both, 
LSMO and LCMO (see the inset of Fig. 2 in Ref.llll). It would 
be highly desirable to have an evaluation which parallels the 
calculations of Millis, Mueller and Shraiman'* for large po- 
larons. For their local modelling they find a temperature de- 
pendent shift of the polaronic resonance for sufficiently strong 
coupling in the undoped case. However, it is not yet conceiv- 
able how such a scheme can be set up for extended excitations 
as the large polarons. 



2. Self-trapped polarons 

The formation of large polarons in the manganites has 
been debated before. Roder, Zang and Bishop'^ studied the 
crossover from a large polaronic state to a quasi-self-trapped 
state for a single hole in the adiabatic limit, applying a vari- 
ational technique for the DE model with the inclusion of one 
effective JT mode. They used a mean-field evaluation for the 
spin degrees of freedom and neglected other electronic corre- 
lation effects. The optical absorption has not been considered 
in this approach. 

Kim, Jung and Noh^ propose that their optical conductivity 
spectra of polycristalline LCMO (Lao.TCao.sMnOa) display 
a sharp Drude peak and a broad MIR absorption band which 
is to be explained as coherent tunnelling and an incoherent 
band of large lattice polarons, respectively. They particularly 
emphasize that, with their electrodynamic analysis, the large 
polaron should be in a strong coupling regime — a polaronic 
state which indeed has not been realized in other physical sys- 
tems because the polaron would usually collapse to a small 
polaron state for strong coupling. However, it is not obvious 
how to identify the large polarons from their optical conduc- 
tivity data. A fit to the large polaron picture by Emin^\ the 
model to which the authors of Ref. !_ refer to, has not been 
presented. The steep low-frequency slope of this type of po- 
laronic absorption cannot be resolved in their data. However, 
since in our thin film optical spectra of LSMO we can clearly 
identify such a threshold behavior, we will now reinvestigate 
Emin's phenomenological approach in some detail. 

The absorption process arises from exciting a carrier from 
the ground state in its (self-trapping) potential well to a free- 
carrier state, invoking the Franck-Condon principle. In this 
intuitive approach, neither the mechanism which is responsi- 
ble for the self-trapping nor the nature of the extended states 
are specified. The optical absorption is calculated by apply- 
ing Fermi's Golden Rule^' which requires the identification of 
the initial and final charge carrier state. Here we slightly de- 
viate from Emin's evaluation by picking a Gaussian localized 
state \G) instead of a hydrogenic ground state as initial state. 
The polaron energy of ]G) is lower and it is better confined 
(see, for example, Ref. 1791). The "Gaussian" result for o'(w) 
will deviate not significantly from that for the hydrogenic 
state. The initial localized electronic wave function is (r|G) = 



(l/i?A/7r)'^^^ exp(— r^/2i?^), and the final free-carrier wave 
function is (r|k) = L~^'^exp(zkr). Emin assumes free- 
carrier states with quadratic dispersion e — k'^ /2m which, 
for the optical absorption, translates into a frequency depen- 
dent wave vector: k{uj) = ^2m{huj — EQ)/h, where Eg/h 
is the threshold frequency for the absorption (Eq is the energy 
difference between the localized ground state and the lowest 
continuum state). The optical conductivity is then: 



a{uj) 



80F< 



1 



n,^^ -(fcHi?)%-M^) 



(9) 



The best fit for the 100 K data is displayed in Fig |8] If 
we assume that m — 3me, consistent with specific heat 
measurements,^'* and Eo/h — 543 cm^^ we find a polaron 
radius from this fit of approximately R/a ~ 2.15/\/3, where 
a ~ 3.88 A. Correspondingly, the diameter of the polaronic 
(ground) state is about 2.5 lattice constants (for m — me, 
R/a ~ 2.15 the diameter is more than 4 lattice constants). 
The maximum height of the resonance fixes the prefactor 
of Eq. (|9ji: the polaron density is estimated for m — Snie 
to a value of iip — 3 x 10^^ cm^"^ (and for m = nie to 
rip — 10^^ cm~'^). The values for R and Up are reasonable in 
this context and support this approach. We point out that the 
hydrogenic ground state, as introduced by Emin and consid- 
ered in our previous work, generates a polaron radius which is 
already unacceptable smalli^^ 

Next we discuss the scaling behavior introduced above (see 
Fig. 13 upper panel). If we rescale a{uj) and oj as 



/ 8y^ Up e^fi 
\ 3 jn Eq 



and w = — , (10) 

Eq 



respectively, and if we also introduce a dimensionless radius 

g = R ^/2mEo/h , (11) 

we identify the rescaled optical conductivity as 
1 



a{w) — — (g ^/w — l) exp(— £1^ (it; — 1)) , (12) 
w 

which only depends on the dimensionless radius g. If g is 
independent of temperature, the observed scaling is also true 
for this approach. For the 100 K data, g — 1.16 ~ ■y/3/2. 
The fits for the other temperatures produce g = 1.10 ... 1.22 
(see TablelllH. a range which is sufficiently narrow to support 
the observed scaling within this approach. 

Even so, is the experimentally estimated value of g ~ 
\/3/2 supported by theoretical considerations? In the 



TABLE III: Self-trapped LP: parameters from Eas.l9ltT2lto fit the 
LSMO MIR spectra at different temperatures for m — Sm^ . 



r(K) 


6 


100 


160 


200 


250 


275 


295 


Eo (cm-') 


526 


543 


580 


650 


777 


876 


1060 


Q 


1.22 


1.17 


1.13 


1.15 


1.14 


1.10 


1.15 


R/a 


1.96 


1.86 


1.73 


1.67 


1.51 


1.37 


1.31 


np (lO^'cm-^) 


2.8 


2.7 


2.6 


2.4 


2.34 


2.26 


2.1 
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Landau-Pekar strong coupling theory with a Gaussian trial 
wave function (see, for example, Ref. 60) the energy 
of the localized polaron is Ep — — (a^/37r) ?ia;LO and 
R = (3\/7r/2/a) i, with the phononic length scale i = 
y/ h / {2muji,o) ■ Correspondingly, the polaron binding en- 
ergy may be reexpressed in terms of the polaron radius: 



p = -(3/2)(a7(2mi?2)). 



If we assumed that \Ep\ = 



Eq holds and if we use the above relation between R and 
g, we would confirm exactly g = \/S/2. Then, since 
a — ^/STTy^Eo/hcJi^o^ we would have for the 6 K data 
1.4 and correspondingly a — 3 ... 4, and 
. 2 with a 



1, 



y/Eo/huji^o 

for the 295 K data ^/E^JKJ^ ~ 1.4 ... 2 with a = 4 ... 6. 
Although this signals the intermediate coupling regime, with 
non-negligible further corrections to the polaron binding en- 
ergy, the estimates seem to be in a reasonable range. 

However, a serious objection is to be raised against this pro- 
cedure: Eq can be identified with \Ep\ only if the self-trapped 
state disappears together with the polarization well in the ex- 
citation process. This should not be the case in the considered 
approach in which we investigate excitations in the adiabatic 
limit and invoke the Franck-Condon principle. In this limit, 
we should take Eq = 3\Ep\ (see Emin's review,^' or the book 
by Alexandrov and Mott,^?5 chapter 1.1). With Eq = 3\Ep\ 
we find g — 3/\/2 which does not fit the experimental datai^ 
Next order corrections to the Landau-Pekar limit do not alle- 
viate this inconsistency. 

Assessment of the strong coupling scenario. The resonance 
in the MIR optical conductivity can be well fitted with Emin's 
phenomenological strong-coupling approach. Again, the fit 
is convincing up to frequencies slightly above the maximum 
whereas the high frequency slope decreases too fast, even 
faster than for the weak coupling scenario. The observed scal- 
ing is approximately reproduced, however with a value for 
the dimensionless radius g which is inconsistent with a strong 
coupling estimate. 

It is mandatory to ask if this phenomenological approach is 
trustworthy in terms of microscopic considerations. Alexan- 
drov, Kabanov and Ray**^ presented numerical results for the 
Holstein model with two and four sites and up to 50 excited 
phonons per mode. In fact, the shape of this microscopically 
calculated polaronic absorption resonance is, for intermediate 
coupling, similar to that of the phenomenological approach: 
a steep slope above the threshold and slower decrease at the 
high frequency side, comparable to the fit curve (Gaussian po- 
laron) in Fig.|S] A direct confirmation of the phenomenolog- 
ical approach was attained for the one-dimensional model. ^^ 
Hence we consider Emin's approach to be reliable for the con- 
sidered adiabatic limit. 

Despite this positive assessment of the theoretical consis- 
tency of the approach we have to question its validity for the 
metallic manganites. The material (LSMO) is conducting and 
the effective mass of the charge carriers is supposed to be 
enhanced but not extremely strong. Consequently, a Fermi 
surface of the carriers should exist above the mobility edge, 
at least for low temperatures (we do not explore any details 
here; for a more comprehensive discussion see the review ar- 
ticle by Edwards^' and the literature cited there, especially 



Refs.|8?, 84). This implies that if such Gaussian localized po- 
larons exist they have to be excited across the mobility edge 
to the Fermi surface — not to a band minimum as assumed in 
the above polaron model. '^^ 

The consequence now is that the phase space for the ex- 
citation of a trapped polaron into "free carrier states" at the 
hole Fermi surface is large, the density of states is finite, and 
one expects a step-like increase of the absorption. The ex- 
perimental shape of the polaronic absorption is, however, not 
consistent with a (possibly smoothed) absorption-edge step. 

Furthermore, one has to account for a mobile charge car- 
rier density of the order of n = 10^^ cm~^ to 10^^ cm^"^. 
This would push the mobility edge^S- to energies well below 
the Fermi edge with the consequence that the frequency of the 
absorption edge is certainly larger than the observed threshold 
— at least for the low temperature measurements. With this 
in mind one should dismiss the scenario that all the polarons 
which contribute to the absorption are self-trapped, that is, lo- 
calized. They may have a large polaronic mass but they are 
mostly delocalized. However, Emin's approach decisively de- 
pends on localized initial states. ^^ For a transition from a band 
state to another free-carrier state one again expects a sharp in- 
crease of the absorption.***^ With these impediments it is un- 
likely that this intuitive approach can be applied directly to 
the observed polaronic absorption in the LSMO films. 



3. Fermi edge singularities 

A generalization of the previous considerations (in 
Sec. IIV C 2> to a finite concentration of mobile carrier states 
is necessary but its implementation is subtle: In the adia- 
batic limit the localized potential well, formed by the lat- 
tice distortion, acts after the excitation of the charge carrier 
to the Fermi edge as a potential scatterer. Consequently, this 
limit represents the x-ray edge problem and the power law 
above the absorption threshold depends on the values of the 
scattering phase shifts.**^ It may diverge (the Mahan exciton 
phenomenon is dominant) or converge (orthogonality catas- 
trophe). Since we have to expect many decay channels for 
this excitation, the singular diverging edge might well be 
smoothed so as to reproduce the shape observed in the mea- 
surements. The decay at the high energy slope is governed by 
ai = 26i{kf)/TT — a where Si{kf) is the phase shift in the 
angular momentum channel I and a is the sum over all phase 
shifts squared (see, for example, Ref. 60). The absorption pro- 
cess involves, according to the standard selection rules, power 
law behavior with both exponents, ai±i. However it is not ev- 
ident, why ai±i should be so strongly temperature dependent 
as to reduce the exponent in the power law to the appropri- 
ate (T-dependent) value. Nevertheless, a might well increase 
with more scattering channels available when approaching a 
disordered paramagnetic state with increasing temperature. 

Above, we argued that we do not expect the polarons, which 
participate in the absorption within the considered frequency 
window, to be localized. With finite polaron dispersion the 
situation changes qualitatively''"'" because the recoil of the 
polaronic potential well has to be accounted for The edge 
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singularity will form a resonance around the frequency of the 
direct transitions from the filled valence band to the Fermi 
edge in the conduction band. This resonance (Mahan exciton) 
will flatten out when the ratio of the valence band mass to the 
conduction band mass increases. It is not clear how to iden- 
tify the polarons with a valence band — or how to quantify 
a "valence band hole" as the appropriate potential well of the 
lattice distortion stripped of its charge carrier. In spite of that, 
the mechanism discussed by Gavoret, Nozieres, Roulet and 
Combescot^" and later on by Ruckenstein and Schmitt-Rinb^i 
is still effective even if the large polarons may not be visu- 
alized as forming a valence band of quasi particles. It is the 
collective behavior associated with the response of the Fermi 
sea in the course of the excitation process that determines the 
optical spectrum. However, we have no information on the 
frequency difference of direct to indirect absorption and we 
have neither any information if the effective mass of the free 
carriers at the Fermi surface increases with respect to that of 
the delocalized polaronic potential well in order to produce 
the observed temperature dependence of the spectra. To sum- 
marize, in a scenario of adiabatic large polarons we have to 
expect these collective effects of the Fermi sea but it is spec- 
ulative if they dominate the shape of the optical spectra in the 
considered frequency window. Furthermore, deviations from 
the adiabatic limit will contribute to a smoothing of the edge 
singularity. 



V. CONCLUSIONS 

We have presented experimental results of the optical con- 
ductivity and dc -resistivity in thin films of LSMO and LCMO 
to analyze polaronic excitations. In both cases the polaron 
formation energy increases with temperature up to Tmi, while 
the spectral weight diminishes continuously. 

Our data indicate different polaronic regimes in both com- 
pounds. For LCMO, the SP picture is well accepted, par- 
ticularly with regard to the high temperature dc-resistivity. 
The optical conductivity displays an asymmetric Gaussian 
line shape close to the maximum of the polaronic resonance 
which also suggests incoherent tunnelling of small polarons. 
Overall, we conclude that for the Ca-doped compound the 
SP model with thermally activated transport is approximately 
valid for both, the ferromagnetic as well as the paramagnetic 
phase. 

In contrast, for LSMO, ln(pfT) deviates from a linear be- 
havior and the optical conductivity has a strong asymmetric 
line shape with a steep slope at the low frequency side. This 
shape is typical for LP absorption processes. -'~^'*" The phe- 
nomenological adiabatic LP model by Emin^' considers ex- 
citations of self-trapped charge carriers to an empty band of 
free-carrier states, with the polaronic potential-well intact in 
the adiabatic limit. Emin's model reproduces the steep low 
frequency slope with a sufficiently wide maximum but cannot 
account for the slow decrease on the high frequency side. We 
challenge that this approach is applicable here because there is 
no hole-band minimum available to which the excited charge 
carriers could scatter. The energy difference from the mobility 



edge, below which the trapped polarons would reside, to the 
hole band minimum would have to be the threshold frequency. 
This frequency Eo/h is (for 6 K) about 526 cm~^, too small 
for the addressed energy difference. Consequently the mod- 
elling has to be modified as to implement LPs with finite mass 
and a finite density of charge carriers (Fermi edge). 

In the weak coupling theory, the generalization of the single 
polaron absorption by GLF'* to the situation of many-particle 
absorption was carried out by Tempere and Devreesei^ The 
optical conductivity for a similar approach, which is rather 
based on the Hubbard model than on a Coulomb gas ansatz, 
is presented in Figs.|8l|9j and II 01 The high frequency decay 
is now better realized through the excitation of particle-hole 
pairs of the Fermi sea. The scaling of the low-frequency slope 
(as observed in Fig.|4] upper panel) is implemented through 
the increase of the number of particle-hole pairs with increas- 
ing energy from the Fermi edge. On the other hand, this model 
cannot account for the shift of the threshold with tempera- 
ture; the model threshold would always be the LO phonon 
frequency which is not supposed to shift from 593 cm^^ for 
6Ktoll83cm-ifor295K. 

The failure of weak coupling theory necessitates to extend 
these considerations to intermediate or strong coupling situa- 
tions, first and foremost in the adiabatic limit. Such an evalua- 
tion is still not available and hardly feasible. As we discussed, 
it involves a collective response of the Fermi sea which leads 
to singularities in the optical conductivity, well known from 
the x-ray edge problem. The recoil of the polaronic lattice de- 
formation, which is taken up in the excitation process of the 
charge carrier to the Fermi edge, will flatten the singularity 
and might be responsible for the observed shape. Although 
the proposed consequences for the shape of the absorption 
spectrum are speculative, the singular response is supposed 
to be present in any model with a sudden creation of a heavy 
or localized scatterer in the presence of a Fermi sea. 

To summarize, the comparison between the measured op- 
tical conductivity spectra and the current theoretical models 
exhibits their inadequacy to explain the spectral weight at the 
high frequency slope of the polaronic MIR resonance and the 
strong shift of the resonance position with temperature. A the- 
oretical solution of the polaron problem which may hold for 
any electron-phonon coupling, including correlation effects, 
is highly desirable. The physical quantities that control the 
temperature dependence for polaron formation, still remain to 
be elucidated. 

Thin films proved to be superior for the analysis of the 
phononic25^ and polaronic spectra in metallic manganites. 
We expect that further investigations based on thin film 
samples will provide considerable additional insight into 
the physics of the metallic state of the CMR manganites. 
We refer to an examination of the doping dependence of 
Lai_j;Sra;Mn03, the polaronic spectra of which may address 
the built up of the Fermi surface. We also propose to explore 
the magnetic field dependence of the thin film spectra in or- 
der to address the role of magnetic correlations in the forma- 
tion of polarons. Finally, Lag. 75(Cao.45Sro. 55)0. 25Mn03 un- 
dergoes an orthorhombic to rhombohedral transition with in- 
creasing temperature"*'' which would invite to examine a pos- 
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sible crossover from SP to LP behavior in the MIR optical 
conductivity. 
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